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Glossaiy 






Null vector 





Null dyadic 


a 


Radius of a spherical subregion 


a, b, di,di 


Arbitrary vectors 


a„ 


(1/2) X maximum linear cross- 




sectional dimension of the sub- 




region V„ 


Atoi 


Dyadic kernel for the MOM 




equations 


Ai,k,Akk 


Self-terms when MOM is used 




with spherical subregions 


b 


Side of a cubical subregion 


B(x) 


Magnetic flux density (weber per 




square meter) 


D{x) 


Electric flux density (coulomb per 




square meter) 


e-mc 


Electric field intensity amplitude 




of an incident plane wave 


E(x) 


Electric field intensity (volt per 




meter) 



<""£«c(x) 


Electric field intensity exciting the 




subregion V„ 


■Ccxc/n 


" ^cxc\Xm ) 


£inc(x) 


Incident electric field intensity 


En, 


Electric field intensity at x„ 


Es^ix) 


Scattered electric field intensity 


/='sca(«0 


Far-zone scattering amplitude in 




the direction «, 


G(x.x') 


Free space dyadic Green's 




function 


Gkm 


= G(xi,x„) 


H(x) 


Magnetic field intensity (ampere 




per meter) 


H^,(x) 


Scattered magnetic field intensity 


i 


= (-1)"^ 


1 


Identity dyadic 


Ax) 


Electric current density (ampere 




per cubic meter) 


Jm 


= Jix„) 


Jconi(x) 


Conduction current density 




(ampere per cubic meter) 
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Jpol(x) 

M 

Pm 

Putn 
'sea 

Qkn, 
t 

u 

u. 

Wine 

Vc« 

Fin, 

Kn 
X 



a 



dm 



8km 

8{x-x') 



Vo 

A 
A 



Polarization current density (ampere 
per cubic meter) 
Free space wavenumber (inverse 
meter) 

Depolarization dyadics for the three- 
dimensional subregion Vm 
Total number of subregions 
Electric dipole moment equivalent to 
the subregion Vn, 
Time-averaged absorbed power 
Time-averaged scattered power 
Time-averaged extinguished power 
Dyadic kernel for the CDM 
equations 

Closed surface of the three-dimen- 
sional subregion Vm 
Time (seconds) 

Polarizability dyadic of the subre- 
gion Vm 

Unit vector denoting direction in the 
far zone 

Unit vector denoting propagation 
direction of an incident plane wave 
Three-dimensional region external to 
the scatterer 

Three-dimensional region occupied 
by the scatterer 
Three-dimensional subregion 
Three-dimensional position vector 
(meter) 

Distinguished point inside the 
subregion Vm 

Xk —Xm 

Mossotti-Clausius polarizability 

scalar of a spherical subregion 

Mossotti-Clausius polarizability 

scalar of the m-th spherical 

subregion 

Kronecker delta 

Dirac delta 

Permittivity of free space 

( = 8.854 X 10"'^ farad per meter) 

Relative permittivity 

Relative permittivity of dielectric 

matter occupying Fim 

Relative permittivity of dielectric 

matter of the subregion Vm 

Intrinsic impedance of free space 

(= IZOvrohm) 

= [(l-/;toa)exp(//:ofl)-l] 

= Wavelength in free space 

( = 2ir/;to) 



Hn Permeability of free space 

( = 4'n-xl0"' henry per meter) 
ft) Circular frequency (radian per 

second) 
n Solid angle 

T Polarizability scalar of a spherical 

subregion 
Tm Polarizability scalar of the m-th 

spherical subregion 
7t) Draine's polarizability scalar of a 

spherical subregion 
7t)B Dungey-Bohren polarizability scalar 

of a spherical subregion 
Vm Volume of the subregion Vm 

(cubic meter) 
i/f(^) A Riccati-Bessel function of 

argument p. 
^(^) A Riccati-Hankel function of 

argument /3. 



1. Introduction 

Electromagnetic scattering problems involving 
complicated geometries are treated numerically 
nowadays. Apart from some low- or high- 
frequency methods [1, 2] and the T-matrix method 
[3], implementation of most numerical techniques 
entails a partitioning of the region occupied by the 
scatterer into may subregions. This is generally true 
whether a differential formulation is used or an 
integro differential formalism. 

The method of moments (MOM) [4-6], as ap- 
plied to an inhomogeneous dielectric scatterer, is 
an approach based on the evaluation of an electric 
field volume integral equation over the region oc- 
cupied by the scatterer. This region is partitioned 
into a number of subregions; the electric field in 
each subregion is represented by a subregional ba- 
sis function; and the volume integral equation is 
converted into a set of simultaneous algebraic 
equations that are solved using standard proce- 
dures. The subregions are generally cubes, 
although the relevant self-terms are usually evalu- 
ated as that of equivoluminal spheres. 

Whereas the MOM is an actual field formalism, 
the coupled dipole method (CDM) is based on the 
concept of an exciting field. The CDM was formu- 
lated intuitively in 1969 by Purcell and Penny- 
packer [7] for dielectric scatterers, although it had 
by that time a rich history spanning many decades 
[8]. Conceived from a microscopic viewpoint, the 
CDM has only a semi-microscopic basis; indeed. 
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the operational basis for applying the CDM to 
boundary value problems is totally macroscopic [9]. 
Both the MOM and the CDM were recently 
extended to bianisotropic scatterers and their 
respective weak and strong forms were shown to be 
equivalent [10]. 

This tutorial paper contains a complete deriva- 
tion of the MOM and the CDM for the scattering 
of time-harmonic electromagnetic waves by an in- 
homogeneous dielectric object possessing an 
isotropic permittivity, the starting point of the exer- 
cise being the monochromatic Maxwell curl equa- 
tions. A central topic of the paper is the 
elucidation of the relationship between the MOM, 
which is widely used in electrical engineering, and 
the CDM, which is motivated by concepts in atomic 
physics [7]. Our treatment is directed towards the 
researcher who is interested in understanding 
these methods but who may not be a specialist in 
electromagnetic field theory. 

A current topic of considerable interest is light 
scattering by agglomerated structures made up of 
individual spheres arranged in a low-density struc- 
ture. Examples of such structures include smoke 
agglomerates formed in fires or internal combus- 
tion engines; various materials produced in com- 
bustion systems including silica and titanium 
dioxide; and, perhaps, interstellar dust. The earli- 
est relevant analyses [11, 12] of light scattering 
from smoke were based on the Rayleigh-Debye 
approximation, in which the field exciting any par- 
ticular primary sphere is taken to be just the field 
that is actually incident on the agglomerate. Such a 
procedure neglects multiple scattering effects, and 
is acceptable for primary spheres with size parame- 
ter ( = radius multiplied by the free space 
wavenumber) less than about 0.2. The typical size 
parameters for the smoke agglomerates mentioned 
above lie in the range 0.1 to about 0.25 at visible 
frequencies. 

Both the CDM [13, 14] and the MOM [15, 16] 
have been applied to compute scattering from 
smoke agglomerates in the recent past, since nei- 
ther technique suffers from the limitations of the 
Rayleigh-Debye approximation. While the MOM 
and the CDM are expected to give the same results 
for infinitesimally small size parameters [17], the 
methods— as ordinarily used— do not yield identi- 
cal results as the primary sphere size increases [18, 
19]. This is because the CDM has been used chiefly 
in what may be called its weak form, while it is the 
strong form of the MOM that is in standard usage 
[10]. The strong form is valid for larger primary 
spheres because the effect of the singularity of the 



free space dyadic Green's function is better esti- 
mated therein than in the weak form. This be- 
comes clearer in the following sections. 

2. Volume Integral Equation 

As is schematically illustrated in Fig. 1, let all 
space be divided into two mutually disjoint regions. 
Kin, and Vm, that are distinguishable from each 
other by the occupancy of matter. The region Kit is 
vacuous; hence, 



D(x) = eoE(xy,x eKcx,, 
B(x) = noH(x);x&Vcr, 



(la) 
(lb) 



The region Km is filled with an isotropic, linear, 
possibly inhomogeneous, dielectric continuum with 
frequency dependent [exp(-/w/)] constitutive 
equations 



D(x) = ene,(jc)£(jc);jceKn., 
B(x) = /^ //(jc); x e Kin, , 



(2a) 
(2b) 



where e,(x) is the complex relative permittivity 
scalar. The square root of e,(x) is the local complex 
refractive index, Imag [weoCr(ji;)] is the local con- 
ductivity, and 0) is the circular frequency. 

There is no requirement that Km be a simply- 
connected convex region. Sharp corners and cusps 
should be absent, it being preferable that the 
boundary surface that separates Kni from Ke« be at 
least once-differentiable everywhere to enable the 
unambiguous prescription of a unit normal at every 
point on it. Furthermore, the maximum linear ex- 
tent of Knt must be bounded so that only the region 
Ken extends out to infinity in all directions. 

The monochromatic Maxwell curi equations, in 
the absence of any externally impressed sources, 
are given in Ke« as 

Vx£(jc)-/w/tt//(jc) = 0;jc e Vrr„ (3a) 

Vx//(jc)+/wco£(4:) = 0;jceKcx,, (3b) 

with denoting the null vector. Similarly, in Kn, we 
have 

S/xEix)-io)tu,H(x) = 0;xe Kn„ (4a) 

VxW(jc) + iwe. €,(x) E(x) = 0;xe Kin,. (4b) 
On rewriting Eq. (4b) as 
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Fig. 1. Schematic of the scattering problem. The unshaded region Vc 
extends to infinity in all directions, while the shaded regions collectively 
constitute V,r,,. 



VxH(x) + ia>eoE(x) + ia>eo [fr(x)-l] £(x) 

= 0;xeH„„ (4c) 

and comparing it with Eq. (3b), we are able to state 
the Maxwell curl equations everywhere as 

'7xEix)-\co^H(x) = 0;xeVi„, + V^, (5a) 

VxHix) + i(oeo Eix)=J(xy, x S Ki„,+ K«. (5b) 

In Eq. (5b), the volume electric current density 
J(x) = for X G Fe„, but 



J(X) = Jpo^x) + J^i{x); X e Fin,, 



(6c) 



which implies that J(x) is not merely a mathemati- 
cally convenient quantity for dielectric scatterers. 
It is enough that we look for the solution of the 
differential equations (5a, b) in terms of only the 
electric field, since the magnetic field everywhere 
may be obtained from Eq. (5a) if the electric field 
is known everywhere: H(x) = VxE(x)/ia>iM). On 
taking the curl of both sides of Eq. (5a), and then 
substituting for Vx/f(j) from Eq. (5b) in the re- 
sult, we get 



J(x)= -ia>€o [e,(x)-l] Eix); x £ Km 



(5c) 



In this fashion a volume electric current density has 
effectively replaced the dielectric matter occupying 
the region Fi„, [20, Sec. 3-11]. 

The volume electric current density /(a:) defined 
by Eq. (5c) is not a fictitious entity in the present 
context. We must remember that, while the region 
Vat is vacuous, the region Ki„, is occupied by dielec- 
tric matter. The monochromatic polarization cur- 
rent density in this dielectric matter is given by 

ypo.(x) = -/&)«« {Real [e,(x)] - 1} E(x); 

X e Ki„,, (6a) 

while the conduction current density is given by 

Jamiix) = CD€o Imag [€,(x)] E(x); x £ Ki„,. (6b) 

It is now easy to see that 



Vx7xE(x)-ko^Eix) = ibifiaJix); 
a:eFi„,-t-Ke«, 



(7) 



where /to = &)(/uoeo)"^ is the free space wavenum- 
ber. We take cognizance of the fact that when J{x) 
is set equal to zero everywhere, Eq. (7) reduces to 
the celebrated vector Helmholtz equation. 

As shown in Appendix A, it follows from Eq. (7) 
that the electric field E{x) is the solution of the 
volume integral equation 

E{x)-EUx)=ia>ii^ fff^^ d'x' 

{G(x,x') -Jix')}; X e Fin,+ Ke„, (8a) 

or, equivalently, 

E(X) -£i„c(x) =ko' iff d'x ' {[€rix ') - 1] 

G(x, X') • Eix')}; X e Ki„, + F„„ (8b) 

on using Eq. (5c). Here 
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G(x , X ') = [I + (l/ito')V V][exp(//to \x -x '!)/ 

4n\x-x'\] (8c) 

is the free space dyadic Green's function and I is 
the identity dyadic. The field E,„c(x) is the solution 
of Eq. (7) when/(jr) = for all x E Vm + V^xx, and is 
the electric field existing in Vinx + Vca if e,(x) = 1 for 
all X £ Fint. The standard radiation conditions are 
satisfied by the right sides of Eqs. (8a,b) as 
/tolrl^oo [21]. 

3. Intermediate Remarks 

Equation (8a) is utilized in setting up the CDM, 
while the integral Eq. (8b) is solved in the MOM, 
and it is this commonality that partially begets the 
algorithmic equivalence of the two techniques. In 
setting up the solutions of Eq. (7) in the form of 
the integral Eqs. (8a,b), we have brought in two 
significant topics that need some rumination right 
now: (i) dyadics and (ii) the free space dyadic 
Green's function. 

Dyadics are as American as apple pie, being the 
brainchildren of Josiah Willard Gibbs. In 1884, 
Gibbs circulated a pamphlet introducing the con- 
cept and nomenclature of dyadics. Mathematics 
books with dyadic notation were written every so 
often earlier in this century, but most mathemati- 
cians appear to have eventually discarded dyadics 
in favor of tensors. In electromagnetics, though, 
dyadic notation has been used with great profit, the 
books by van Bladel [22], Fedorov [23] and Chen 
[24] being immensely popular. A short exposition 
on dyadics was brought out by Lindell [25] in 1981, 
and is much recommended to the interested 
reader. 

A dyadic serves as a linear mapping from one 
vector to another vector; thus, a dyadic D is a map- 
ping from a to ft given by ft = D • a . This property 
leads to the idea of a dyad that is composed of two 
vectors, i.e., Dn=did2. It follows that 
Di2'a =di{di'a) and a -012 = (a •di)d2 are vectors, 
and Di2Xa =di(d2Xa) and a xDi2=(a y.d\)di are 
dyads, and the oft-used appellation bivectors for 
dyads appears justified [26]. Because a dyadic can 
be written as the sum of dyads, the general repre- 
sentation of a dyadic is the sum D = .5*m Dkmdi4m 
[27]. 

The identity dyadic I is such that D • I = D = l • D, 
and the null dyadic is such that D • = = 'D. 
The simplest antisymmetric dyadic is « x I, where « 
is any vector of unit length. Even vector differential 
operators can be thought of as dyadics; thus, the 
curl operator is written as V x I, and the divergence 



operator as V • I, in dyadic notation. It is not possi- 
ble for us to go through all the wonderful proper- 
ties of dyadics than can be exploited in theoretical 
electromagnetics research, so we refer the inter- 
ested reader to the compendiums in the books by 
Chen [24], van Bladel [22], and Varadan et al. [28]. 
The main feature of computational significance is 
that, as dyads are bivectors, all dyadics used in this 
paper can be thought of as 3 x3 matrices. 

The second issue concerns the singularity of the 
dyadic Green's function G(jr, x') used in Eqs. 
(8a,b) and defined in Eq. (8c): the factor exp (//to 
|jir-jr'l)/47rljr-jf'l goes to infinity as ljr-jr'l-*0. 
When evaluating G(jr, x'), the second derivatives 
therefore have to be carefully handled. As shown 
by van Bladel [29], the classical procedure leads 
one to write 

///^d'x'G(x,x')-ft(x')=-(l/3y)ft(x)-h 

P.V. USy d'x' {G(x, x')'b{x')); xEV, (9) 

where P.V. stands for "principal value." When 
computing the P.V. integral on the right side of Eq. 
(9), an infinitesimally small spherical region cen- 
tered about the point j:'=x is excluded from the 
domain of integration. Computing the P.V. integral 
is not problematic because [30] 

Gi{x,x') = [{\-uxUx) + {ilk,\X\)(\+ilh\X\) 

{S-'iuaix)] {&xp{ikM\)l4iv\X\};\X\^Q, (10) 

is not singular when the source point x'and the 
field point x do not coincide; here X=x -x' and 
ux=XI\X\. Yaghjian [31] has modified the right 
side of (9) to an expression wherein the excluded 
infinitesimal region need not be spherical. 

We will use an alternative approach to evaluate 
the integral (9), as shown in Appendix 2. In this 
approach [6, 10, 32], the region of integration is 
split into two regions: one region includes the point 
x'=x in its interior, while the second one is the 
remainder. Use is made of Gauss' theorem and the 
Green's function for Poisson's equation to deter- 
mine this integral. This procedure is attractive as it 
places very little restrictions on the shape and the 
size of the region surrounding the point x' =jr. 



4. The Method of Moments 

The method of moments is a general mathemati- 
cal procedure for transforming a linear operator 
equation into a set of simultaneous algebraic 
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equations, and the interested reader is referred to 
[4], Chap. 1, for a simple introduction. We are, 
however, confining ourselves here to the applica- 
tion of MOM to electromagnetic scattering prob- 
lems. 

Although the MOM has grown increasingly so- 
phisticated in the last decade [6, 33], a simple ver- 
sion [5] suffices for the easy conversion of the 
integral equation (8b) into a set of simultaneous 
algebraic equations. We begin by partitioning the 
scatterer region Fi„t into simply connected subre- 
gions Vm, Ism <M, each bounded by a closed sur- 
face 5m, so that Eq. (8b) can be transformed to 

Eix)=E,„,(x) + k„' £m SSL d'x' {Mx')-1] 



G(x,x')-E(x% X S X„V^ + Vc«. 



(11) 



The main features of this partitioning scheme are 
as follows: 

(i) Each subregion Vm is modeled as being 
homogeneous such that 

£r(x) = £r,m;xE Vm. (12) 

(ii) The maximum linear cross-sectional extent 
2a„ of V„ is such that flm/A<0.1 and 
«mler,m"^l/A <0.1; that is, the dimension a^ is 
no more than a tenth of the wavelength A in 
the exterior region Vex, as well as of the wave- 
length in the subregion Vm. 

(iii) The surface 5m that bounds the subregion V^ 
is sufficiently smooth so as to be at least once 
differentiable, which enables the unambigu- 
ous prescription of a normal at any point on 

Satisfaction of these three conditions, in prac- 
tice, means that the union SmVm of the subregions 
is only approximately coincident with the scatterer 
region Vm. Condition (i) can lead to an artificial 
material discontinuity across the interface of two 
adjacent subregions, therefore the simple MOM 
algorithm given in this section works best when . 
adjacent subregions do not differ widely in their 
permittivities. Condition (ii) ensures that the spa- 
tial variations of the electromagnetic field inside 
each subregion are small enough so that each sub- 
region can be thought of as a dipole scatterer [34], 
though different authors use somewhat different 
upper bounds on the subregional size [35, 36]. Con- 
dition (iii) is mostly ignored by MOM users, their 
usual practice being to use cubical subregions. 
Thus, the adequacy and the accuracy of the 
MOM— and the CDM— results depend in large 
measure on the adequacy of the partitioning 
scheme. 



Two ancillary aspects of the partitioning scheme 
need to be stated for the simple MOM algorithm. 
First, the incident field must have slow enough spa- 
tial variations that it may be considered almost spa- 
tially constant over any subregion. Second, not only 
the permittivity but also the actual field is assumed 
spatially constant in each subregion. Together, 
these two assumptions constitute a long wavelength 
approximation [37], whose consequences are ex- 
ploited in the MOM as well as in the CDM. 

Let the volume of Vm be denoted by Vm = JJJy 
d^x, and let Xm denote a distinguished point (such 
as the centroid) lying inside Vm. On setting jf =jfi in 
Eq. (11) we get the relation 

E(x,)-ko'[e,,,-l]fSS^^d'x' 

{G(X, , X') ■E(x')}=Ei„ciXk)+ko^ Zn.m>k W„ - 1] 



//L d?x- {G{x,, X') -Eix')}; 1 <k ^M. 



(13a) 



Next, on using the approximation E{x)=E{xm) for 
all xEV„ and carrying out the integrations over all 
subregions, we get 

£t-(ea-l)[VM* -U]-Ek = 

Eir,c(Xk) + ICo Sm.mfk {Vm [fr.m " 1] 



G(Xk,x„,)'E„}; \<k<M, 



(13b) 



where Em = E(x„,). In obtaining Eq. (13b) from Eq. 
(13a), we have accomplished the integration on the 
subregion V„,, m T^k, very simply by evaluating the 
specific integrand at Xm and multiplying it by the 
volumetric capacity Vm. The integral ///^, die' 
G{Xk,x')'E{x'), on the other hand, has been esti- 
mated using Eq. (A2-5a). 

With the assistance of Eqs. (5c) and (12), we 
transform Eq. (13b) into 3A/ simultaneous alge- 
braic equations compactly stated in vector dyadic 
notation as 

Ei„^{xk) = lm^(\x-.A() [A*m -Em]; Is k <.M , (14a) 

where the 3M unknowns are the cartesian compo- 
nents of the fields £m, ISm <M . The dyadic kernel 
used in Eq. (14a) is given as 

Atm = I 5i„, -^ (1 - Er,* ) 5im [^O^M* -\.k] + 

ko'vm (1 - €r.m ) (] - 5*m ) G*m , (Mb) 
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where 8km is the Kronecker delta and the definition 

Gkm = G(Xk, X„) = [(I -A*^*m/IAto, I') + 
{ilk^km\y {\-3Xkn,Xlanl\Xlan\'-)\ 



[e\^(ika\Xkm\)IAir\XkJ]; m T^k, 



(15) 



follows from Eq. (A2-6a) with Xkm =Xk —Xm. In this 
straightforward and simple version of the MOM, 
wherein the electric field has been assumed to be 
piecewise constant over the scatterer region Vi„,, 
Eq. (14a) has to be solved for the three cartesian 
components of all Em for specified E\nc{x). 

Bearing in mind that all of our dyadics can be 
easily thought as 3 x 3 matrices, one can solve Eq. 
(14a) using a variety of matrix manipulation tech- 
niques. The Gauss elimination procedure [38] is 
simple but places a heavy demand on computer 
memory. Much less computationally intensive is 
the conjugate gradient method [39], which is now 
being heavily employed for MOM calculations [40]. 

Once all Em have been thus calculated, the scat- 
tered electric field in Fe,, can be computed as the 
sum 

EUx)=E{x)-EUx) = 

ko' Xme{\,l....M){Vm [Cr.m - 1] G{x,Xm) • E„]\ X 6 Vm, 

(16) 

which expression follows from Eq. (8b). Now, from 
Eq. (A2-6a), we have the asymptote [41] 

Tend,„^.,„,^.G(x,x„) = [\-{x-x„){x-Xm)l 

\x-\n\^] [exp(/folx -xJ)I^tt\x -x„\\, (17a) 

correct to order 1/lx -x„\. Since allxm are generally 
distributed around the origin, the further reduction 

Tend,„,.,„,_ G(x,x„.)= [l-xc/UfP] 

{exp(/A:olj:-x„l)/4T7lxl}, (17b) 

is reasonable. Finally, making the Fraunhofer ap- 
proximation [42] 
exp(/A:olx-x„l)sexp[//:olJfl] exp[i/coX •x„/Lcl] (17c) 

on the left side of Eq. (17b) yields 



Tendi^^.,^,^. G(x, Xm)s [l-xc/brl^] exp(//:ol*l) 

exp(-/A;ox-x„/lxl)/4T7|x|]. (17d) 

Consequently, the scattered electric field of Eq. 
(16) has the asymptotic representation 

Tendkin-.- E^{xu^= [e.\p{iknK)lx] F«^(us), (18) 

in the direction «s in the far zone, with the far-zone 
scattering amplitude defmed by 

F«;a(«0=(^oV47r) [|-«s«s]. 

^m^[,i..M) {exp{-ikMs-x„) v„ [e,.„ - 1] E„}, (19) 
which shows quite clearly that «s •Fsca(«s) =0. 



5. The Coupled Dipole Method 

The heart of the MOM is Eq. (13b) which in- 
volves the electric field Ek that is actually present at 
Xk . However, in the CDM we have to consider the 
electric field that excites the subregion Vk, each 
subregion being explicitly modeled as an electric 
dipole moment. In order to obtain an expression 
for this exciting field, we return to Eq. (8a), and 
rewrite it as 

E(x)-ia>^ //4 d'x' {Qix,x')-J{x')} = 

Ei„cix) + i(OfM, lm.m*k llly^ d'x' {G(x,x') J (x')}; 

X6X.Kn+F.«, (20) 

after partitioning the scatterer region precisely in 
the same manner as was done for the MOM solu- 
tion. We know that £i„c(x) is the field in the ab- 
sence of the scatterer and that the quantity 

i<^\y^ //4„ d'x' {G(x, x') •y(x')}; x^V„,m^k, 

is an electric field whose source clearly lies in the 
subregion Vm, m^k. Thus, the whole of the right 
side of Eq. (20) is a composite electric field with 
multiple sources, but none of the sources lies in Vk. 
This composite electric field 

<*'£«c(x)=£inc(x) + /wA10 lm.m*k U ly^ ^^X' 

{G{x,x-)-J{x')\,x&Vk, (21a) 

can be thought of as the field that excites Vk. 
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Since the left side of Eq. (20) must be equal to 
its right side, it follows that 

<*^^cxc(x)=£(x)-/a)pfl //^^ d'x' {G(x,x')'J (x')}; 

xen. (21b) 

On setting x =Xk in Eq. (21b), and completing the 
volume integration therein using Eq. (A2-5a), we 
get the expression 

E^,t = Ek -iafJo [Mt-(Vkti^)Lk] Jk =Eac.k; 

i£k<M; (22) 

here,/i =J(xk) and E„ck =***£exc(*t). But, by virtue 
of Eq. (5c) we have 

Jk =i(oeo (1 - ea ) Ek ; 1 s A: <M, (23a) 

which means that Eq. (22) reduces to 

Aw • Ek =E,^k ■,\<k<,M, (23b) 

whence, 

Ek =Akk-' 'E^,k\ \^k^M, (23c) 

Jk =/fc»6o (1 - er,i) Am "' • £.xc.i ; 1 sA: <A/. (23d) 

Consistent with the long-wavelength approach, in 
the CDM we think of the equivalent electric dipole 
moment [9, 10] 



Pk =(i/(o)vkJk 
located atx*; thus, 

Pk=tk 'Ecxc,k, 

after using Eq. (21b), where 

U=vt 63(e,,i-l) Aw~' 



(24a) 



(24b) 



(24c) 



is the polarizability dyadic of the dielectric region 
Vk. 

We return to Eq. (21a) at this stage, set x =Xk 
therein, and do the volume integrations; the result 
is the vector dyadic relation 

Ecxe,k =E\oc(Xk)+i(OIXo Xn.m^k {Vm Gfo,, -Jm); 

\<,k<M, (25) 



in accord with the long-wavelength approximation. 
With the aid of Eq, (24a), we rewrite Eq. (25) as 

£e«c,t =£inc(r*) + w^/io Xm.m^k {Gfe, '/>„} , (26a) 
and the use of Eq. (24b) yields 

£«ci=£:i„c(Xi)-f-fc)^Pfl Xm.m*k{Gkm-\m}-E„c^. (26b) 

Equation (26b) constitutes the core of the CDM 
and, for numerical work, is best rewritten as the set 
of the 3A/ algebraic scalar equations 

E\ac{Xk) = Xma{\,X...M) {Q)m' EeK,m\, \^k^M, 

(27a) 

for the cartesian components of £«K.m, where the 
dyadic kernel 

Qhn = \8k„,-o}^yi^Gkm-U{\-8km) . (27b) 

The CDM algorithm is structurally just the same as 
the MOM algorithm, as is easily demonstrated by a 
comparison of Eqs. (27a,b) with Eqs. (14a,b). 
Hence, computational techniques for solving Eq. 
(27a) are the same as for solving Eq. (14a). 

Once all £„c,mhave been obtained as the solution 
of Eq. (27a), one can use Eq. (24b) to find all pm 
and ascertain the scattered field as the sum 

E^{x)=E{x)-E-„^(x) = 



(0 



'iMi Xmc {1.2... M) {G(r, Xm) 'pn,}; xEVcxt , (28) 



which relation follows from Eq. (8a), Equation (18) 
still applies for the scattered electric field in the far 
zone, but the far-zone scattering amplitude is now 
given by 

Xnc{i.2...M) {exp(-ikou, 'x„) p„}; (29a) 

equivalently, 

^«:a(«s) = ((O^/Mi/^Tt) [I -M5H5] ' 
-SnG{l,2.. , M\ {eXp( -ikoUs •Xn,)ln- £cK.m}. (29b) 

Parenthetically, we remark here that the 
Rayleigh-Debye approximation for scattering by 
smoke agglomerates [11, 12] may be obtained by 
replacing £«c,m in Eq. (29b) by E;„,(x„ ). 
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6. Scattering and Absorption 

The presence of a scatterer has two observable 
consequences insofar as the conservation of elec- 
tromagnetic energy is concerned. A portion of the 
energy radiated by the source of £inc(jr) is scattered 
in all directions by the scatterer. Another portion is 
absorbed by the scatterer and converted into other 
forms such as heat. All calculations pertaining to 
the scattered and the absorbed powers at a given 
circular frequency w are made in terms of power 
averaged over the time period 1ttI(o. 

Because Fxa{ui) is of the form [I -KsKs] 'b,'n fol- 
lows that Us'Fsca(us) = 0, as has been mentioned 
earlier. By virtue of Eqs. (3a) and (18), we have 

Tend*or-.„ Hsaixus)= 

iVrjo) [eTip(ikox)/x] UsXFsc=,ius), (30) 

where 7jo = (/W6o)"^ is the intrinsic impedance of 
free space. These observations imply that the scat- 
tered field in the far zone is Transverse Electro- 
Magnetic (TEM) in character. 

The time-averaged scattered power per unit 
solid angle in the direction Us is defined as 

dP«a(Ms)/d/3(K,) = (l/2) 

Tendio^^, Real [x^u,'{E^,(xu,)xH^*(xus)}],(31a) 

with dn(Ui) = s\nei dft d0s, as is customary in 
spherical coordinates, and the asterisk denoting 
the complex conjugate. After putting Eqs. (18) and 
(30) in Eq. (31a), we obtain 

dPsca(Ms)/d/2(Ms) = (l/27Jo)Fsca(K,) " /^s«*(k,); (31b) 

consequently, we are able to ascertain the total 
time-averaged scattered power as the integral 

/'sca = (l/2Tio) I d^ I dftsin^s 

F^.(us)-F^*(Us). (31c) 

Unless the scatterer material is intrinsically loss- 
less, there is absorption of electromagnetic energy 
in Fint. Because the scatterer is simply a dielectric 
object here, the time-averaged power absorbed in 
ViM may be computed as the volume integral 



from Poynting's theorem for monochromatic fields 
[24, Chap. 2]. After using Eq. (2a) on the right side 
of this relation, we get the result 

/'ab, = Real[(/w/2) 

//4„ d'x {e^,*ix)Eix) • £*(x)}] . (32b) 

Finally, after using Eq. (12) as well as the long- 
wavelength approximation made heretofore, we are 
able to reduce this expression to the sum 

Pab,=Real [(/&)/2) //f . d'j:{€o€r,,„* £„ ■£„*}] = 

' inl 

- (w/2) Imag [X, u„{€o€r.m *E„'E„. *}] . (33) 

Because E„ • E„* is purely real, no power absorp- 
tion will take place in V„ if e,^ is purely real; in- 
deed, Pabs^O, provided Imag [6r(jc)]=0 for all 

X e Vin,. 

Insofar as the MOM is concerned, the solution 
{E„; m = 1, 2, ..., M} of Eq. (14a) may be directly 
substituted into Eq. (33) for the computation of 
Pabs. The calculation of Pabs in the CDM is only 
slightly more complicated: the exciting fields 
{EcK^; m = \, 2, ..., M} obtained by solving Eq. 
(27a) have to be substituted into Eq. (23c) to get 
{E„;m = 1, 2, ..., M}. Thus, an useful expression for 
CDM computations is 

Pabs = - (W /2) Imag [X. Vm{£66r.m * 

'm ' i^cxc^ ' 'm 'iicxc^ /J , W / 

where 

T„ =[! + (! -6,„) (kMn.-l^)]-' . (34b) 

The total time-averaged power extinguished by the 
scatterer region is the sum 



1 cxt — ^sca "t" J abs . 



(35) 



Pabs=Real [(/w /2) fff^, d'x £(jr) •Z)*(jr)] 



(32a) 



Quite often, one is interested in the extinction of 
the plane wave 

Ei„c(x) =einc exp[ikou\„c -x] , (36a) 

Hir^ix) = (I/tjo) Wine X Ci« exp[i/:(iMj„c -x] , (36b) 

where e^c carries the units of volts per meter and 
Mine is a dimensionless unit vector such that 
einc*«inc = 0. In this case, the total time-averaged 
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power extinguished by the presence of matter in 
V\M can be estimated using the forward scattering 
amplitude as [22, Chap. 8]: 

Pcxt = (2Tr/w/io) Imag [Cmc* • Fsca(Minc)] . (37) 

Substitution of Eqs. (29a) and (36a) in this relation 
gives us 

Pc^ = {co/2) Imag [X„{Ei^*(x„) -p^}] , (38a) 

whence, 

Pc« = {(0/2) Imag [Xn{iBi„c*(x„) • t, • £c.c^}] (38b) 

for use with the CDM. The MOM analog of Eq. 
(38b) is easily obtained by substituting Eq. (23b) 
therein; ergo, 

Pcxt = ((oeo/2) Imag [X.{v„,(er,m -1) 

E;„,*(x„,) ■£„}], (38c) 

for use with the MOM. 

7. Spherical and Cubical Subregions 

Although spheroidal and ellipsoidal subregions 
have been used [43, 44], it is commonplace in liter- 
ature to have cubical or spherical subregions. CDM 
users are more comfortable with spherical subre- 
gions, while MOM users are fond of cubical ones. 
Cubes and spheres have the same dyadic L, and in 
many MOM codes the M dyadic of a cube is esti- 
mated as that of an equi-voluminal sphere; see 
Appendix C. Without any particular loss of gener- 
ality therefore, we take the subregions to be spheri- 
cal and of identical radii in the remainder of this 
paper. 

Let the subregion Vi, be the sphere of radius a 
with its center at x*. As a result, the volume 
Vk =(4^^/3)fl^ the dyadic Lt = (1/3)1, and the dyadic 
Mk=(2/3ko^)[(\-ihra) exp(*oa)-l]I, as shown in 
Appendix C. Consequently, the MOM dyadic ker- 
nel given in Eq. (14b) reduces to 

A«=l[l + (l/3)(l-ea){2(l-/A:,^) 

exp(//toa)-3}], (39a) 

Afon =(477/3) a^ ko' (1 -er.m)G*„; m ^k . (39b) 

In the same manner, the CDM dyadic kernel of Eq. 
(27b) simplifies to 



Q« = l, (40a) 

Qkm = -(o^fMi i4v/3)a^ eo(er/7. - 1) 

[l-^(l/3)(l-e^,m){2(l-/A:ofl) 
exp(ik(fl ) - 3}] - ' Gfan ; m ?=* . (40b) 

An analysis of the self-terms is crucial to the 
understanding of the MOM and the CDM. From 
Eq. (39a) it follows that the MOM self-term can be 
broken up as 



where 



and 



Akk={Auc+Auc}\ , 



/l«=(er,t + 2)/3. 



(41a) 



(41b) 



A =(2/3)(l -ea)[(l -ikifi) exp(/A:,^)-l] . (41c) 

Both Avt and Akk should be called self-terms; in- 
stead, it seems only Akk has been accorded that 
honor in some MOM papers, e.g., [16-19, 45]. 

The CDM self-term is somewhat obscure, being 
hidden as the polarizability dyadic tt of Eq. (24b). 
In the present instance, the polarizability dyadic 
reduces to a scalar because the subregion is spheri- 
cal; hence, 

\k = \rk, (42a) 



where 



and 



rk=akl{\+AkklAkk), (42b) 



0* = 4Trfl ' eo(er.t - l)/(er,t + 2) (42c) 

is the Mossotti-Clausius polarizability [46, 47] of 
the electrically small dielectric sphere. 

Let A^ofl < < 1 in Eq. (41c) and Akk be evaluated 
correct to order h?a^. Then, 

Tk =akl{\ -ko-(a -' + 2ikJ3)ak/4Treo} , (43) 

and we observe that the (2/feV3)oit/4Treo term in the 
denominator of the right side of Eq. (43) is the 
radiative reaction term used by Draine [48] in his 
CDM formulation. More often, Akk is evaluated 
correct only to the first order in k(fl, leading to 
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Tk = cn^ [7], and thereby giving rise to the semi- 
microscopic flavor of this numerical approach [9]. 

8. Strong and Weak Forms 

The various approximations that can be made for 
the self-terms lead us to the strong and the weak 
forms of the MOM and the CDM [10]. In the 
strong forms (S-CDM & S-MOM), Eq. (A2-5a) is 
used to estimate the self-term (m =A:) in Eqs. (13a) 
and (21b). In the Weak forms (W-CDM & W- 
MOM), Eq. (A2-7) is used in place of Eq. (A2-5a) 
for this estimation. For isotropic dielectric scatter- 
ers, the W-CDM is exemplified by Purcell and Pen- 
nypacker [7] using spherical subregions, and the 
S-CDM has become recently available [17] for the 
same problem. The S-MOM has been used for 
isotropic dielectric scatterers for many years [6, 
33], but the idea of W-MOM is of more recent 
provenance [17]. 

The W-MOM corresponds exactly to the 
W-CDM, while the S-MOM corresponds exactly to 
the S-CDM. When all subregional volumetric ca- 
pacities are very small, the S-MOM/S-CDM effec- 
tively transmutes into the W-MOMAV-CDM. 
Generally stated, therefore, it can be surmised that 
the scatterer region Fim must be broken up into a 
larger number of subregions when the W-MOM/ 
W-CDM is used than if the S-MOM/S-CDM is 
used. Comparison of S-MOM results with the 
W-CDM results, with identical partitioning of the 
scatterer region, [e.g., 16, 18, 19], in some instances 
may be like comparing apples with oranges. The 
difference between S-CDM/S-MOM and W-CDM/ 
W-MOM is solely due to the inclusion or the exclu- 
sion of the dyadic M* in the expression for A«, 
therefore computational time is marginally in- 
creased and the memory requirement negligibly, 
when one shifts from W-CDM/W-MOM to 
S-CDM/S-MOM. 

9. Spherical Subregions and Finite-Size 
Effects 

An assessment of the MOM and the CDM for 
isotropic dielectric scatterers with spherical subre- 
gions is now in order. To facilitate such a compari- 
son, we reiterate that 

Att = {H- (1 - Cr.* )[(2/3)(l - fArofl ) 

exp(;^(^)-l]}l, [S-MOM] (44a) 

A,, = {(ea + 2)/3}l ; [W-MOM] (44b) 

and correspondingly. 



n = (47r/3)fl ' co(ea - 1)/{1 + (1 - e^ ) [(2/3) 

(1 - ik^ ) expiiki^ ) - 1 ]} , [S-CDM] (45a) 

TicS Ava^ eo{e,,k - l)/(6a +2) = Qt . 

[W-CDM] (45b) 

We can evaluate the strong forms vis-a-vis the 
weak forms through the Taylor expansion of the 
ratio r/a about Cr = 1, there being no need for us to 
continue the subscript k in the remainder of this 
paper. Using Eqs. (41b), (41c) and (42b), we get 
the expansion 



T/a = l+2vl[(er-l)/3] + 
2A(,2A- l)[(e,-l)/3]2 + 



(46a) 



where A = [(1 -ikua) exp(JA:ofl)-l] has been used 
for convenience; and it becomes clear that the 
weak forms are poor approximations of the strong 
forms when the relative permittivity is significantly 
different from unity. 

Likewise, the Maclaurin expansion of r/a with 
respect to the size parameter koa. 



T/a = l+[(e,-l)/(e, + 2)](A:,^)2 + 
+ (2z73) [(e,- l)/(e,+2)](M)'+ .... 



(46b) 



indicates that the weak forms become increasingly 
inadequate as the size parameter of the subregion 
increases. 

Equations (46a) and (46b) tell us that strong 
forms incorporate finite-size effects meaningfully, 
while the weak forms do so trivially. This conclu- 
sion is reinforced by Figs. 2-4 that show plots of 
rfa versus the normalized radius Arofl of an electri- 
cally small sphere for 6,= 1.5, 2.5 and 4.0. In these 
figures we have ensured that the maximum value of 
^ofllcr"^ is about 0.5. We observe -not surprisingly 
—that T is complex while a is purely real for these 
input parameters. Furthermore, in all three figures 
the ratio Ir/al lies between 1.02 and 1.03 when 
^ufll6r"^s0.5; while T/a = l for k<fi =0, of course. 
Similar conclusions can be drawn from Fig. 5, 
wherein the calculations of r/a have been reported 
for a lossy dielectric sphere with 6, = 3.75 -t- 0.25/. 
These figures highlight the understanding on finite- 
size effects drawn from Eqs. (46a,b), and it follows 
that coarser partitioning of Ki„, may be acceptable 
when S-CDM/S-MOM is used than when W-CDM/ 
W-MOM is used. These thoughts can be validated 
by careful examination of the graphs recently 
published by Ku [19]. 
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Fig. 2. Functions [T/or-l] (solid lines) and [Tbii/a-l] (dashed 
lines) as functions of the size parameter lufi of a spherical sub- 
region whose relative permittivity £,= 1.5; the Mossotti-Clausius 
polarizability a = 47ra-' eu{e,-\)/(e, + 2). 



Fig. 4. Same as Fig. 2, but €, = 4.0. 
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Rg. 3: Same as Fig. 2, but €, = 2.5. 



ng. 5. Same as Fig. 2, but €, = 3.75 + 0.25i". 
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The strong forms discussed above are not the 
only way of introducing finite-size effects, there be- 
ing at least two more CDM algorithms available for 
that purpose. As alluded to earlier, Draine [48] 
replaced the Mossotti-Clausius polarizability by 

7b = a/{l-(i/67reo)io'a}. [D-CDM] (47) 

However, td is fully contained in t, as we have seen 
in Eq. (43); and td approaches a as the size 
parameter ktfl becomes vanishingly small. 

The second source for the incorporation of 
finite-size effects has been the Lorentz-Mie-Debye 
analysis for scattering by a dielectric sphere [46]. 
The electromagnetic field scattered by any object 
can be described in terms of a multipole expansion 
[49]. With inspiration from Doyle [50], Dungey and 
Bohren [51] introduced finite-size effects in the 
CDM by using the lowest order Mie coefficient— 
corresponding to the electric dipole term of the 
mutipole expansion — for the polarizability in place 
of the Mossotti-Clausius polarizability a. Thus, 



Tt>B=(i6lT€o/ko^) 



[DB-CDM] (48) 

was used in [47], with i/f(i3) = j8"' sin(j8) -cos(j8), 
^i/^()8) = di/^/d)8, f()8)=-(l+i/)8) exp(;)8), and 
^f()8) = df/d^. 

Using Eq. (48), it can be shown that tdb ap- 
proaches a as the normalized radius k^ becomes 
vanishingly small. This is also borne out in Figs. 
2-5, wherein the computed values of Tte/a are 
compared with those of r/a as functions of the nor- 
malized radius k(ff for er = 1.5, 2.5, 4.0 and 
3.75 + 0.25/. The general behavior of tub is the 
same as of t: (i) tob is complex even for purely real 
er, and (ii) the ratio iTtje/al lies between 0.99 and 
1.01 when ^ofl ler"^l s 0.5 in these figures. A more 
remarkable observation is that II - Tbe/al S II - rla\; 
in other words, tdv is closer to the Mossotti- 
Clausius polarizability a than t is. 

Draine [48] and Dungey and Bohren [51] 
concluded from their numerical investigations that 
D-CDM and DB-CDM, respectively, generally 
provide scattering results superior to those from 
W-CDM. This does not come as a surprise because 
the self-terms in W-CDM (or W-MOM) are 



estimated with the least accuracy. On the other 
hand, although it is difficult to provide general 
enough conclusions for the adequacy of either 
D-CDM or DB-CDM vis-a-vis that of the S-CDM/ 
S-MOM, it is safe to state that any claims of superi- 
ority—based purely on the estimation of some 
gross parameter, such as the total scattering cross- 
section— are debatable. Indeed, the only good way 
of deciding on the superiority of S-CDM, D-CDM 
or DB-CDM is by making calculations for the 
specific problem under consideration: The scat- 
terer region should be subdivided into different 
numbers of subregions, and computed data from 
the various algorithms for different partitioning 
schemes compared for the property of interest [52- 
54]. 

The refractive index of soot at visible frequencies 
has been measured by a variety of experimental 
techniques, and has been found to be dependent 
on the source materials [55]. We conclude with 
calculations of r/a and TDe/a for £, = 2.40 + 2.38/, 
corresponding to a complex refractive index of 
1.7 + 0.7/, in Fig. 6. It is not difficult to gather from 
this figure that conclusions drawn in the previous 
paragraphs of this section apply to light scattering 
by soot agglomerates. 



0.03 




Fig. 6. Same as Fig. 2, but t, = 2.40 + 2.38/. 
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10. Appendix A. The Volume Integral 
Equation for the Electric Field 

The free space Green's function G(x,x') satisfies 
the dyadic differential equation 

VxVxG(x,x')-ko'G(x,x') = 

\S(x-x'), (AM) 

where S(x-x') is the Dirac delta function. Since 
J(x') is not a function of x, it follows that 
Vx[G(x,x')'J(x')] = [VxG(x,x')]-J(x'); 
hence, 

VxVx[G(x,x')-J(x')]- 

ko^ [G(x,x')-J(x')] = l-J(x') S(x-x') . (Al-2) 

We multiply both sides of this equation by iw/M) and 
integrate over all x' to get 

(vxvxi-kon)-fff^,^^d'x' 

[i(oiJLoGix,x')-J(x')] = 

'■'"MO ///,,,, d'x'J(x') S(x -X') , (Al-3) 
whence, 

(VxVxl-Wl)-ni,.,^.d'x' 

[/wMo G(x,x') -/(r')] =icotu>J(x) . (Al-4) 



that depends on the forcing function, J(x) being the 
forcing function here, and the complementary 
function that holds when the forcing function is 
identically zero everywhere. The right side of Eq. 
(Al-6) is the particular solution of Eq. (7), which 
has the complementary function Emc(x) satisfying 

^x VxEi^(x)-ko%„c(x)=0; 

X 6Fi„, + Kc«. (Al-7) 

Hence, the complete solution of Eq. (7) is 

E(x)^Ei„c(x)+i(oiJUi fffy,^ d^x' 

[G(x,x') -Jix ')]; X e Ki„, + Kx, , (Al-8) 
which can be easily re-expressed as Eq. (8a). 

11. Appendix B. The Self-Integral 

The integral under consideration is given as 

a (xo) = //L d^ ' G(xo, x')-b(x'), (A2-1) 



where Kis the region bounded by the surface 5, as 
shown in Fig. 7, G(x, x') is given by Eq. (8c), and 
the point xo lies inside V. Because G(x,x') is of the 
order l/lr -x'\^ for Ix -x'\= 0, and becomes singu- 
lar at x=x', this integral has to be carefully 
Now by comparing Eqs. (7) and (Al-4), it is not too treated. 



arduous to convince ourselves that 

£(x)=/6,Mo ///,.,,. d^' [G(x,x')'J(x')]; 

OsUrl<oo. (Al-5) 

Remembering that for the present purposes 7(x) is 
null everywhere except in Vm, we next obtain 

E(x)=i(oixo fff^,^ d'x' [G(x, x') -/(x')]; 

X eKi„, + K«. (Al-6) 

The solution of every linear differential equation 
can be divided into two parts: ih& particular solution 



As shown by Fikioris [56], we can transform Eq. 
(A2-1) to 

fl(xo) = ///^.^_^d'x'[G(xo,x')-6(x')] + 



//4^ d'x' [G(xo, X') •*(x')-G,(xo, X') •ft(xo)] 



-{\lko^)fLdh'u„' 
[(x ' -Xo)/47rlx ' -xol'] • b (xo) , (A2-2) 



where 

G,(x,x') = (l/Jto') ^^[l/47rlx-x'l] (A2-3) 
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is an auxiliary dyadic related to the Green's func- 
tion for Poisson's equation [10]; Vo is an exclusion- 
ary region bounded by the surface So, as shown in 
Fig. 7; u„ is the unit normal to So, pointing away 
from Vo, at the point x' £5o. The exclusionary re- 
gion Vo should be small but not necessarily in- 
finitesimal, and it must be wholly contained within 
V. Moreover, there is no requirement that So be a 
miniature copy of S. 

The use of the long-wavelength approximation 
b(x)sb(xo) for all x £ F reduces Eq. (A2-2) to 



a(xo)=SSS^_^^d'x' [G(xo,x')'b(xo)] + 



SSSyd'x'{[G(xo,x')-Gsixo,x')]-b{xo)} 



■■Ko 



-(l/ko')JLd'x'u„'[(x'~Xo)/ 



'So 



4'!T\x'-Xo\']-b(XQ) 



(A2-4) 



Since V is electrically small, but Vo need not be 
infinitesimal, we set V=Vo and 5 =So to have 



a (xo) = [M - O/kd") ^ • b (xo) , {A2-5a) 



where 



M = //I d'x ' [G{xo, X ') - a (xo, X ')] , {A2-5b) 



L = ;;^ d^ ' «„ ' [{x ' -xo)/47rlx ' -xol'] . {A2-5c) 

The evaluation of M can be accomplished nu- 
merically for regions of different shapes using the 
coordinate-free (i.e., non-indicia!) representations 
[24, Chap. 9] 

G(x, X ') = [(I -KAMA") + {ilko\X\){\ + ilh\X\) 



(1-3 uaix)\ [exp {iko lx-x'l)/47rlx-x'l]; 



\X\i^Q, 



(A2-6a) 



a(x,x') = (i/^oi:^l)' (1-3 UxUx)\l4'u\X\- 

\X\^Q, (A2-6b) 




v-v. 



VoE 



Fig. 7. For the evaluation of the integral in Eq. (A2-1) when Xa E Fand 
x'BV. 
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v/heTeX=x-x' and ux=X/\X\. The evaluation of 
the depolarization dyadic L is equally easy on a dig- 
ital computer, and analytical expressions for L are 
available for the quite general ellipsoidal shapes 
[57]. 

If Vis extremely small in electrical size, the sim- 
pliflcation 



a{x^) = {l/ko^)L'b{xo) 



(A2-7) 



is permissible and leads to the weak forms of the 
MOM and the CDM. 



12. Appendix C: The Self-Dyadics for 
Spherical and Cubical Regions 

Let the region V in Appendix B be a sphere with 
radius a and center at xq. Without loss of generality 
due to the translational invariance of the right 
sides of Eqs. (A2-6a,b), we set Xq=0. Furthermore, 
we set x' = t' («i sinfl' cos<^'+«2 sinfl' sin0'+«3 
COS0'), where lx'l=r', (u],U2,U3) is the triad of 
cartesian unit vectors, while (r', 0', <t>') represent 
the corresponding spherical coordinate system lo- 
cated at the center of the sphere. Hence, 



!SSyd'x'{-}^\ dr'r'^ I dfl'sinfl' f 'd0'{-}, 
Jo Jo Jo 



ffs^hc' {'} =a^ I de' sine' { d<^'{-}. 
Jo Jo 



Now, 



'U}filx = 



I de' sine' [ d<^i 
Jo Jo 

f de'sine'f d0'{(«, sin0'cos0' + 
Jo Jo 

«2 sine' sin<^' + «3 cos0')(«i sinfl' cos<^'+«2 
sinfl' sin<^' + «3 cos0')} = 

dO' sine' {(«i«i-l-«2«2) sin^0-l-2«3«3 cos^0'} = 
Jo 

(47r/3) («i«i+«2U2+«3«3) = (47r/3)l, (A3-1) 



and 



f dfl'sinfl' f d<^' l=47rl (A3-2) 

Jo Jo 



likewise. It follows that 



M 



Upherc= dr'r'^ dfl'sinfl' d<^' 
Jo Jo Jo 

{[G(0,x')-a(0,x')]} 

= f dr' r" [(21/3) exp(fA:or')/r'] 
Jo 

= I {mU) {{\-ik^) exp(fA:^ ) - 1] . (A3-3) 

In order to evaluate L for a spherical region, we 
note that u„' =x'l\x'\ where x' E. S; ergo, \x'\=a, 
and 

Ll,phere=fl^ f dB' sinB' \ d<t>' 
Jo Jo 

{«„'x'/477a^} = (l/47r) 

J de'sine'J d<^' {(«i sinO'cosf^'-t- 
Jo Jo 

«2 sinfl' sin<^'+«3 cos0')(«i sinfl' cos<^'-l- 
«2 sine' sin<^'+«3 cosfl')} 
= (l/47r) (47r/3) I = (1/3)1 , (A3-4) 



after using Eq. (A3-1). 

If the region V in Appendix B were to be a cube 
of side b, we would get 



Llc.be = (1/3) I, 



(A3-5) 



the same as that for a sphere. No expression exists 
for the dyadic M for a cube in closed form; however, 
it is commonplace in MOM practice to use the 
value for an equivoluminal sphere. Thus, 

Micobe = I (2/3A:o') {[1-ikJ) {S/AttY^] 

exp[«A:,*(3/4)"']-l}, (A3-6) 



714 



Volume 98, Number 6, November-December 1993 

Journal of Research of the National Institute of Standards and Technology 



because [(3/477)"^ b] is the radius of a sphere hav- 
ing the same volume as the cube of side b. 



Acknowledgments 

One of the authors (George W. Mulholland) was 
supported in part by the Pennsylvania State 
University Particulate Materials Center while on 
sabbatical leave at Pennsylvania State University. 
The authors thank Dr. Egon Marx (NIST), Prof. 
Dike Ezekoye (University of Texas at Austin), and 
Prof. W. S. Weiglhofer (University of Glasgow) for 
their constructive comments and criticisms. Thanks 
are also due to the external reviewers whose com- 
ments helped improve this tutorial paper. 



13. References 

[I] H. Massoudi, C. H. Dumey, and C. C. Johnson, Long- 
wavelength electromagnetic power-absorption in ellip- 
soidal models of man and animals, IEEE Trans. 
Microwave Theory Tech. 25, 47-52 (1977). 

[2] G. I. Rowlandson, and P. W. Barber, Absorption of 
higher-frequency RF energy by biological models: Calcu- 
lations based on geometrical optics. Radio Sci. 14 (6S), 
43-50 (1977). 

[3] V. V. Varadan, A. Lakhtakia, and V. K. Varadan, Com- 
ments on recent criticism of the T-matrix method, J. 
Acoust. Soc. Amer. 84, 2280-2284 (1988). 

[4] R. F. Harrington, Field Computation by Moment Meth- 
ods, McGraw-Hill Book Publishing Company, New York 
(1968). 

[5] D. E. Livesay, and K. M. Chen, Electromagnetic fields 
induced inside arbitrarily shaped dielectric bodies, IEEE 
Trans. Microwave Theory Tech. 22, 1273-1280 (1974). 

[6] J. J. H. Wang, Generalized Moment Methods in Electro- 
magnetics, John Wiley and Sons, New York (1991). 

[7] E. M. Purcell, and C. R. Pennypacker, Scattering and ab- 
sorption of light by nonspherical dielectric grains, Astro- 
phys. J. 186, 705-714 (1973). 

[8] A. Lakhtakia, General theory of the Pureell-Pennypackcr 
scattering approach, and its extension to bianisotropic 
scatterers, Astrophys. J. 394, 494-^99 (1992). 
[9] A. Lakhtakia, Macroscopic theory of the coupled dipolc 
approximation method. Opt. Commun. 79, 1-5 (1990). 

[10] A. Lakhtakia, Strong and weak forms of the method of 
moments and the coupled dipolc method for scattering of 
time-harmonie electromagnetic fields, Int. J. Mod. Phys. 
C 3, 583-603 (1992); errata: 4, 721-722 (1993). 

[II] M. V. Beny, and I. C. Percival, Optics of fractal clusters 
such as smoke. Optica Acta 33, 577-591(1986). 

[12] R. D. Mountain, and G. W. Mulholland, Light scattering 
from simulated smoke agglomerates, Langmuir 4, 1321- 
1326 (1988). 

[13] J. Nelson, Test of a mean field theory for the optics of 
fractal clusters, J. Mod. Opt, 36, 1031-1057 (1989). 

[14] S. B. Singham, and C. F. Bohren, Scattering of unpolar- 
ized and polarized light by particle aggregates of different 
size and fractal dimension, Langmuir 9, 1431-1435 (1993). 



[15] M. F. Iskander, H. Y. Chen, and J. E. Penner, Optical 
branching and absorption by branched chains of aerosol, 
Appl. Opt. 28, 3083-3091 (1989). 

[16] J. C. Ku, and K. H. Shim, A comparison of solutions for 
light scattering and absorption by agglomerated or arbi- 
trarily-shaped particles, J. Quant. Spectrosc. Radiat. 
Transfer 47, 201-220 (1992). 

[17] A. Lakhtakia, Towards classifying elementary microstruc- 
tures in thin films by their scattering responses, Int. J. 
Infrared Millim. Waves 13, 869-886 (1992); errata: 14, 
663 (1993). 

[18] J. I. Hage, and J. M. Greenberg, A model for the optical 
property of porous grains, Astrophys. J. 361, 251-259 
(1990). 

[19] J. C. Ku, Comparisons of coupled-dipole solutions and 
dipole refractive indices for light scattering and absorp- 
tion by arbitrarily shaped or agglomerated particles, J. 
Opt. Soc. Amer. A 10, 336-342 (1993). 

[20] R. F. Harrington, Time-Harmonic Electromagnetic 
Fields, McGraw-Hill Book Publishing Company, New 
York (1961). 

[21] D. S. Jones, Low frequency electromagnetic radiation, J. 
Inst. Maths. Applies. 23, 421-447 (1979). 

[22] J. van Bladcl, Electromagnetic Fields, Hemisphere, 
Washington, DC (1985). 

[23] F. I. Fedorov, Theory of Gyrotropy (in Russian), Nauka i 
Tcknika, Minsk (1976). 

[24] H. C. Chen, Theory of Electromagnetic Waves, McGraw- 
Hill Book Publishing Company, NcwYork (1983). 

[25] I. V. Lindell, Elements of Dyadic Algebra and Its Appli- 
cation in Electromagnetics, Radio Laboratory Report S 
126, Helsinki University of Technology (1981). 

[26] P. Boulanger, and M. Hayes, Electromagnetic plane 
waves in anisotropic media: An approach using bivectors, 
Philos. Trans. Roy. Soc. Lond. A 330, 335-393 (1990). 

[27] L. Brand, Vector and Tensor Analysis, John Wiley and 
Sons, New York (1947). 

[28] V. V. Varadan, A. Lakhtakia, and V. K. Varadan, eds., 
Field Representations and Introduction to Scattering, 
North-Holland, Amsterdam (1991). 

[29] J. van Bladcl, Some remarks on Green's dyadic for 
infinite space, IRE Trans. Antennas Propagat. 9, 563-566 
(1961). 

[30] W. Weiglhofer, Delta-function identities and electromag- 
netic field singularities, Amer. J. Phys. 57, 455-456 
(1989). 

[31] A. D. Yaghjian, Electric dyadic Green's functions in the 
source region, Proc. IEEE 68, 248-263 (1980). 

[32] J. J. H. Wang, A unified and consistent view on the singu- 
larities of the electric dyadic Green's function in the 
source region, IEEE Trans. Antennas Propagat. 30, 463- 
468 (1982). 

[33] E. K. Miller, L. Mcdgycsi-Mitschang, and E. H. Newman, 
eds.. Computational Electromagnetics: Frequency- 
Domain Method of Moments, IEEE Press, New York 
(1991). 

[34] D. S. Jones, Numerical methods for antenna problems, 
Proc. Inst. Elcc. Eng. 121, 573-582 (1974). 

[35] K. M. Chen, and B. S. Guru, Inlcmal EM field and 
absorbed power density in human torsos induced by 
1-500-MHz EM waves, IEEE Trans. Microwave Theory 
Tech. 25 746-756 (1977). 

[36] M. J. Hagmann, O. P. Gandhi, and C H. Dumcy, Upper 
bound for cell size for moment method solutions, IEEE 
Trans. Microwave Theory Tech. 25, 831-832 (1977). 



715 



Volume 98, Number 6, November-December 1993 

Journal of Research of the National Institute of Standards and Technology 



[37] C. Athanasiadis, Low-frequency electromagnetic scatter- 
ing by a multi-layered scatterer, Quart. J. Mech. Appl. 
Math. 44, 55-67 (1991). 

[38] B. Camahan, H. A. Luther, and J. O. Wilkes, Applied 
Numerical Methods, John Wiley and Sons, New York 
(1969). 

[39] G. Strang, Introduction to Applied Mathematics, 
Wellesley-Cambridge Press, Wellesley, MA (1986) Chap. 
5. 

[40] T. K. Sarkar, ed.. Application of Conjugate Gradient 
Method to Electromagnetics and Signal Analysis, El- 
sevier, New York (1991). 

[41] We define g{x) = Tendt^„ f{x) as a function that repli- 
cates the asymptotic behavior of f{x) as j;-»», "Tend" 
standing in for "tendency;" Lim^^. [&(*)//(•>;)] = 1- 

[42] J. W. Goodman, Introduction to Fourier Optics, 
McGraw-Hill Book Publishing Company, New York 
(1968) Chap. 3. 

[43] S. B. Singham, Intrinsic optical activity in light scattering 
from an arbitrary particle, Chem. Phys. Lett. 130, 139-144 
(1986). 

[44] V. K. Varadan, A. Lakhtakia, and V. V. Varadan, Scatter- 
ing by beaded helices: Anisotropy and chirality, J. Wave 
Mater. Interact. 2, 153-160 (1987). 

[45] G. H. Gocdecke, and S. G. O'Brien, Scattering by irregu- 
lar inhomogeneous particles via the digitized Green's 
function algorithm, Appl. Opt. 27, 2431-2438 (1988). 

[46] C. F. Bohren, and D. R. Huffman, Absorption and Scat- 
tering of Light by Small Particles, John Wiley and Sons, 
New York (1983). 

[47] Chronological propriety demands that Mossotti's name be 
placed before that of Clausius; see: M. T. Prinkcy, A. 
Lakhtakia, and B. Shanker, On the extended Maxwell 
Gamett and the extended Bruggeman approaches for 
dielectric indiclectric composites, Optik (1993; in press). 

[48] B. T. Draine, The discrete dipolc approximation and its 
application to interstellar graphite grains, Astrophys. J. 
33, 848-S72 (1988). 

[49] C. H. Papas, Theory of Electromagnetic Wave Propaga- 
tion, Dover Publications, New York (1988) Chap. 4. 

[50] W. T. Doyle, Optical properties of a suspension of metal 
spheres, Phys. Rev. B 39, 9852-9858 (1989). 

[51] C. E. Dungey, and C. F. Bohren, Light scattering by non- 
spherical particles: a refinement to the coupled-dipolc 
method, J. Opt. Soc. Am. A 8, 81-87 (1991). 

[52] H. Massoudi, C. H. Dumcy, and M. F. Iskander, Limita- 
tions of the cubical block mode! of man in calculating 
SAR distributions, IEEE Trans. Microwave Theory Tceh. 
32, 746-752 (1984). 

[53] V. V. Varadan, A. Lakhtakia, and V. K. Varadan, Scatter- 
ing by three-dimensional anisotropic scatterers, IEEE 
Trans. Antennas Propagat. 37, 800-802 (1989). 

[54] R. A. West, Optical properties of aggregate particles 
whose outer diameter is comparable to the wavelength, 
Appl. Opt. 30, 5316-5323 (1991). 

[55] I. Colbeck, E. J. Hardman, and R. M. Harrison, Optical 
and dynamical properties of fractal clusters of carbona- 
ceous smoke, J. Aerosol Sci. 20, 765-774 (1989). 

[56] J. G. Fikioris, Electromagnetic field inside a current- 
carrying region, J. Math. Phys. 6, 1617-1620 (1965). 



[57] A. Lakhtakia, Polarizability dyadics of small chiral ellip- 
soids, Chem. Phys. Lett. 174, 583-586 (1990). 



About the authors: Akhlesh Lakhtakia is Associate 
Professor of Engineering Science and Mechanics at the 
Pennsylvania State University. George W. Mulholland 
is Research Chemist and Group Leader, Smoke 
Dynamics Research Group at NIST. The National 
Institute of Standards and Technology is an agency of 
the Technology Administration, U.S. Department of 
Commerce. 



716 



